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ABSTRACT 

For explaining the broadband spectral shape of knots in the M87 jet from 
radio through optical to X-ray, we propose a modified synchrotron model that 
considers the integrated effect of particle injection from different acceleration 
sources in the thin acceleration region. This results in two break frequencies at 
two sides of which the spectral index of knots in the M87 jet changes. We discuss 
the possible implications of these results for the physical properties in the M87 
jet. The observed flux of the knots in the M87 jet from radio to X-ray can be 
satisfactorily explained by the model, and the predicted spectra from ultraviolet 
to X-ray could be further tested by future observations. The model implies that 
the knots D, E, F, A, B, and CI are unlikely to be the candidate for the TeV 
emission recently detected in M87. 

Subject headings: galaxies: active — galaxies: individual (M87) — galaxies: jets 
— radiation mechanisms: non-thermal 



1. Introduction 

The giant radio galaxy M87, situated nearly at the center of the Virgo cluster at a 
distance of 16 Mpc (e.g. Macri et al. 1999), has been studied extensively from the radio 
(e.g., Owen et al. 1989 and references therein; Biretta et al. 1995; Sparks et al. 1996), 
optical (e.g. Perlman et al. 2001a, hereafter POl; Biretta et al. 1991), UV (Waters & Zepf 
2005, hearafter WZ05), to X-ray (e.g., Marshall et al. 2002, hereafter M02; Wilson & Yang 
2002). The predicted TeV gamma-ray (Bai & Lee 2001) and the possible position where 
the TeV emission from the M87 jet (Georganopoulos et al. 2005; Cheung et al. 2007) have 
also been conflrmed at a high confidence level by the HESS observations (Aharonian et al. 
2003 & 2006). In particular, Perlman & Wilson (2005, hereafter PW05) analyzed diagnostics 
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and physical interpretation of the X-ray emissions from the M87 jet in detail. Many other 
observations have also been done, such as photometric surveys of the jet, polarimetry maps 
of the jet (Perlman et al. 1999). 

In the M87 jet, the Optical- X-ray spectral index ctox (-5',^ oc increases from 1.2 — 1.4 
at 2" - 7" from the nucleus (knots D and E) to 1.7 - 1.9 at 15" - 18" from the nucleus 
(knots B and CI). The Chandra data confirm steep X-ray spectra of the knots in the M87 
jet (Ox > ttr) and arc thus consistent with a synchrotron origin for the X-ray jet emission 
(M02; Wilson & Yang 2002; PW05). So far, there have been three standard theoretical 
synchrotron models. The KP model (Kardashcv 1962, hereafter K62; Pacholczyk 1970) 
assumes that the source of the emission is a single burst of energetic electrons with an 
isotropic pitch-angle distribution and thus no scatters. Because of the likely scattering of 
relativistic particles by hydromagnetic waves (e.g., Wentzel 1977), the KP model is physically 
unreasonable, and therefore will not be mentioned further in this paper. The JP model ( Jaffe 
& Perola 1973) assumes the same initial conditions as those of the KP model, but allows 
scattering in the pitch-angle distribution so that it can maintain an isotropic distribution all 
the time. The resulting spectrum is an essentially exponential rollover above the synchrotron 
loss break frequency. The continuous injection (CI) models (K62; Heavens & Meisenheimer 
1987, hereafter HM87) assume that a power law distribution of relativistic particles is being 
continuously added to the emitting region, but the CI model of HM87 further takes the 
advective transport of the accelerated electrons downstream into account. The CI model of 
K62 has the similar spectral shape to the one of HM87. 

But these synchrotron models also have some problems. As shown in WZ05, considering 
the X-ray flux, the CI model is only apphcable for the inner knots (knots D and E), and 
explains the UV turnover (Fig 1), but for other outer knots (such as knots F and A) this 
model systematically over-predicts the observed UV turnover. And PW05 also found that 
such a model cannot explain the spectral index at X-ray. Without considering the X-ray 
data, the predicted X-ray flux by the CI model is higher than the observed one, but the 
exponential high-energy rollover of JP model underpredicts the X-ray flux by many orders 
of magnitude and the slope at X-ray is much larger than the observed. The two standard 
models can't fit the X-ray flux and the spectral index at X-ray under the single index of the 
electron energy spectrum at injection and single emission process. We summarize the criteria 
for fitting the spectral energy distributions (SEDs) of the knots in the M87 jet discussed in 
WZ05 and PW05: firstly, the first break frequency should be under the UV turnover (Fig 1), 
especially, for knots D, A, and B; secondly, a steeper X-ray spectral index than the optical is 
needed, so there may be a second break frequency between UV data and X-ray data; finally, 
the best fitting model should explain the fiux and index of the X-ray data as well as the ones 
of the radio, optical and UV data. 
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In § 2, we describe in detail our modified synchrotron model for knots in the M87 jet. 
In § 3, we present and discuss the fitting results of this model to the wide-band spectra of 
the knots D, E, F, A, B, and CI. A summary is given in § 4. 



2. The Model 

We suggest a synchrotron model which is a modified CI model of K62, to explain the 
observed SEDs of the knots in the M87 jet. The kinetic equation of the relativistic electrons 
in the CI model (K62) is : 

?^ = ,1.(E^„)^,E-', (1) 

where (3 — hB\, 6 is a constant, represents the component of the magnetic field perpen- 
dicular to the velocity of the particle. The synchrotron energy losses arc dE/dt = —^E"^. 
The injection of a constant spectrum qE'^ in the CI model of K62 actually implies that 
the acceleration region is the same as the main emission region in the knots. However, this 
assumption may not be true in general, especially for shocks acceleration mechanism which 
may be the dominant mechanism of the particle acceleration in jets. Based on the observa- 
tional evidence of shocks in the M87 jet (Perlman et al. 1999; PW05; Harris & Krawczynski 
2006) , a more physical scenario would be that the acceleration region and the main emission 
region are not strictly co-spatial in the M87 jet, and the accelerated injection electrons may 
be advected downstream (this process is similar to the CI model of HM87). The electrons 
advection throughout the jet and diffusion throughout the jet's cross-section with the de- 
crease of the synchrotron lifetimes from low energies to high energies may result in spatially 
stratified emission regions along the jet (M02) and the consistent narrowing of the jet from 
radio to optical to X-ray (PW05). For the acceleration region, a more complex but physical 
scenario may be that there are many acceleration sources in a compact acceleration region 
which is unresolved by the telescope beam and each source has a power law distribution 
of injection population of relativistic electrons. We assume that all the sources accelerate 
electrons under the similar conditions (e.g., similar magnetic fields etc.) with the same mech- 
anism which results in a same spectral index p, similar electron number density, the same 
maximum energy and fiuid velocity for the injection populations of relativistic electrons into 
the main emission region. Before injecting into the main emission region, the populations of 
relativistic electrons in the acceleration region will experience a brief synchrotron loss pro- 
cess. But the observed knot emission is dominated by the emission from the main emission 
region. If the injection electrons into the main emission region come from a very compact 
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region, not resolved by the observation, we need to integrate all the contributions of popu- 
lations from different locations of the acceleration region by considering synchrotron losses. 
We assume that r is the distance to the main emission region from a source of the acceler- 
ation region. At low energies, the populations from different locations of the acceleration 
region have the same spectral shape, so the integrated injection spectrum of the electrons 
over the acceleration region still scales as E~p. At higher energies, the electrons are advected 
a distance r to the main emission region before losing a significant fraction of their energy to 
synchrotron emission. As the loss time-scale scales as E'^, this distance scales as E~^. Thus 
the integrated injection spectrum at high energies would be proportional to E~^'^^\ There is 
a break energy E^ between the low energies and the high energies in the integrated injection 
spectrum. The integrated injection spectrum exhibits a cut-off energy Ec- By adding up the 
contributions from sources in the acceleration region, the integrated injection spectrum as 
the source into the main emission region from the acceleration region can be written as, 

qiE-P, EkE^; 

q,E-^+'\ E,>E>Eb; (2) 
0, E>Ec, 



where, q^E'^^^^ corresponds to the component of the injection spectrum at high energies. 
The spectrum near the break energy in fact is not a strict power law, and the equivalent 
spectral index should be intermediate. But here we would simply consider a power law form 
for the spectrum. The real distribution of the electrons and the magnetic field in the M87 
jet is very complex (M02; Perlman et al. 1999; Perlman et al. 2001b), but we simply assume 
that the distribution of the electrons is isotropic and the large-scale equipartition magnetic 
field in the main emission region is approximately symmetrical. The kinetic equation is: 



dt ^dE^ ' 



qiE-P, E<Eb; 
q2E-^+^\ E>Eb. 



(3) 



Then, we have 



N{E,e,t) 



gig-(P+l) 

I3{p-1) ' 



p-ii 



^ < jl ! 
E,>E>f^ 

E> Eb. 



(4) 



For isotropic distributions, we can derive the flux expression of the modified synchrotron 
model 
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U > Ub2 ■ 



(5) 



r^Bi = c^t-\ VB2 = C3(3^El C3 = 3.4 x 10^5-3, 



where t (in yr) is synchrotron hfetime, ubi and 1/^2 (in Hz) are the first and second break 
frequencies, respectively. 



Now, we apply the above modified synchrotron model to the observed knot emission in 
the M87 jet. These include knots D, E, F, A, B, and CI (as showed in Figure 1 of Harris & 
Krawczynski 2006). 

The data we used are listed in Table 1. These include the published radio to UV data 
(POl), UV data at 1.8 x 10^^ Hz (WZ05) and X-ray data (PW05). There is another reported 
X-ray measurement (M02) made 12 days before PW05 observation. M02 did not use the 
same regions to integrate the fluxes for various components as did POl. PW05, however, did 
use the same regions as POl, so we use PW05 data points in our model fitting. There was no 
X-ray measurement for the knot CI. So we estimated its flux density using the measurements 
for the knot C and the ratio of knots CI to C from the optical observation at 1.0 x 10^^ Hz 
by assuming this ratio is the same at both optical (1.0 x 10^^ Hz) and X-ray. There were 
no error bars in the X-ray data by PW05. The error bars for the PW05 X-ray data listed 
in Table 1 were estimated by assuming the same relative precision (the uncertainty in terms 
of a fraction of the value of the result) in both PW05 and M02 data. We use the weighted 
least square method to fit our model (Equation 5) to the wide-band data in each knot, in 
this progress a power law form is chosen near the break frequencies because we only need 
to concern the trend of the break frequencies. There are two peak frequencies in our model, 
but we don't know where they are when we fit our model to the broadband data in each 
knot. So we first arbitrarily divide broadband data into three groups to perform the least 
square method, we could calculate the corresponding reduce chi square {xD by changing the 
division. All the possible combinations are considered before we obtain the best fit with a 
minimal xl among them. These best- fit parameters for each knot, are listed in Table 2. The 
SEDs for the knots in the M87 jet and the best fits for our model are plotted in Fig 1. 

Our model can well fit the radio, optical, UV, and X-ray data for these knots with three 
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segments except knot E which has two segments implying that the second break frequency of 
knot E has exceeded the X-ray band. It satisfies the aforementioned three constraints for the 
SED of the knots of M87 jet (as shown in Fig 1 and Table 2). The predicted spectra shape 
of knots from UV to X-ray (from —p/2 to —{p + l)/2) could be tested by future telescope 
in the band. The equivalent spectral index of the high energy electrons in the acceleration 
region may not decay from p to p + 1, but is likely to be between them (i.e. the high energy 
electrons may be near the break energy) . This results in an equivalent spectral index of the 
photons possibly between p/2 and {p-\- 1)/2 (e.g. knot D). Our model is consistent with the 
spectral indexes of the knots E, F, A, B, and CI within their error bars. 

The first break synchrotron frequencies {vbi) of all the knots (D, E, F, A, B, and CI) 
in the M87 jet are under 10^^ Hz (see Table 2). According to the unified scheme of AGNs 
and equation (1) in Bai & Lee (2001), the Compton components from these knots will peak 
at 0.01 ~ 1 GeV, and their fiux densities at TeV would be undetectable (based on the 
synchrotron sclf-Compton model and the equation (4) in Bai & Lee 2001). This implies 
that the knots D, E, F, A, B, and CI arc unlikely to be the candidate for TcV emission in 
the M87 detected by the HESS observations (Aharonian et al. 2003 & 2006). The possible 
candidate for the TeV emission in M87 is related to the innnermost region of M87 like HST-1 
or core itself (Stawarz et al. 2006; Aharonian et al. 2006; Cheung et al. 2007). HST-1, the 
knot closest to the core of M87, has been shown to have a very dynamic fight curve and 
flaring (Perlman et al. 2003; Harris et al. 2003; Stawarz et al. 2006; Harris et al. 2006). Its 
mechanism may be related to the core and thus more complicated. 

From Table 2, as a whole the second break frequency vb2 decreases along the jet. This 
may imply that the synchrotron loss in the acceleration region will increase along the jet. 

We find that the value of the particle spectral index p is about 2.36 on average, which 
agrees well with the latest expectations from both diffusive shock acceleration theory (2.0- 
2.5, Kirk & Dendy 2001) and acceleration by relativistic shocks (2.23 in the ultrarelativistic 
fimit. Kirk 2002). Particle acceleration at shocks (e.g., Blandford & Ostriker 1978) through 
the first-order Fermi process is generally believed to occur in jets. 

PW05 proposed a possible phenomenological model to modify the CI model. They 
suggest that the filling factor of the volume within which particles are accelerated declines 
with increasing energy at X-ray energy along the jet (not in the direction perpendicular 
to the jet), but they can't explain the running mechanism of filling factor. To explain the 
curved X-ray spectra of BL Lac objects, Perlman et al. (2005) consider an episodic particle 
acceleration model which assumes only a time-variable particle acceleration. This results in 
a logarithmic curvature rather than a sudden break and could be related to the broadband 
spectral shape too. Fleishman (2006) suggests a very different model which explicitly takes 
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into account the effect of the small-scale random magnetic field, probably present in the 
M87 jet. But the energy densities contained in small-scale and large-scale magnetic fields 
may be incomparable, we think that the electrons in the large-scale magnetic field could also 
give rise to emission of the knots in our model, especially at high frequencies. The idea of a 
secondary population of the relativistic electrons that have different spectral index from the 
first population is discussed by Jester et al. (2005), which is partially similar to our model. 
But we consider a lot of populations of relativistic electrons that have the same spectral 
index in the acceleration region and discuss the detailed process that may be responsible for 
the knots in the M87 jet. 



4. Conclusion 

We propose a modified CI model that considers a decay of spectral index of injection 
electrons possibly due to the sum of the injection spectrum from different acceleration sources 
with synchrotron losses in the thin acceleration region, so there are two break frequencies at 
two sides of which the spectral index changes for the spectra of knots in the M87 jet. We 
consider that the emission of the knot may be still emitted by the relativistic electrons in the 
large-scale magnetic field at high frequencies as well as the low frequencies, but not by the 
small-scale random magnetic field (e.g., Fleishmann 2005). Our model gives a satisfactory 
fit to the SEDs of knots in the M87 jet. Based on our analysis, the knots D, E, F, A, B, and 
CI are unlikely to be responsible for the TeV emission detected in M87. The fitting results 
from our model imply that the particles in M87 are accelerated by shocks, and as a whole 
the second break frequencies of knots decrease down the jet. We also predict the spectra of 
knots from UV to X-ray, which could be tested by future observations in the band. 
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Table 1. Flux Densities of Knots in the M87 Jet 



Flux Density 



Frequency 

V (Hz) 


D 


E 


F 


A 


B 


CI 


Observing 
Date 


Reference 


VLA (mJy): 


















1.50 X 10^° 


161.54 ±1.92 


48.05 ± 0.81 


144.90 ± 1.86 


1218.00 ±12.00 


808.40 ± 8.30 


544.70 ± 5.60 


94 Feb 04 


POl 


HST (^Jy): 


















1.45 X 10^'' 


280 ± 11 


71.4 ±3.0 


262 ± 11 


2633 ± 105 


1739 ± 70 


971 ± 39 


98 Apr 04 


POl 


1.87 X 10^* 


224 ±9 


67.2 ±2.0 


271 ±11 


2344 ± 90 


1489 ± 60 


797 ± 32 


98 Feb 26 


POl 


2.66 X 10" 


168 ±7 


42.9 ±2.3 


147 ±6 


1829 ± 73 


1104 ±44 


606 ± 24 


98 Feb 26 


POl 


3.75 X 10" 


150 ±4 


42.3 ± 1.3 


158 ±5 


1363 ± 34 


811 ±20 


415 ± 10 


98 Feb 25 


POl 


4.50 X 10" 


117±3 


33.4 ±0.9 


123 ±3 


1086 ± 27 


623 ± 16 


310.7 ±7.8 


98 Feb 25 


POl 


6.58 X 10" 


96.2 ±2.5 


28.0 ±0.8 


98.8 ±2.6 


904 ± 23 


504 ± 13 


241.0 ±6.1 


98 Feb 25 


POl 


1.00 X 10^^ 


59.5 ±1.6 


16.2 ±0.6 


62.7 ± 1.7 


586 ± 15 


306.8 ±7.8 


135.9 ±3.5 


98 Feb 25 


POl 


1.80 X lO^'^ 


26.7 ±0.9 


10.0 ±0.9 


34.5 ± 1.7 


275.2 ± 12.7 


147.4 ±7.9 


80.8 ±12.9 


01 Feb 23 


WZ05 


Chandra (njy): 


















2.42 X 10^^ 


51.5 ±4.2 


32.2 ±6.5 


20.1 ±5.2 


156 ±8.8 


30.3 ±5.5 


14.6^ ±5.2 


00 Jul 29 


PW05'' 



"^The X-ray flux density of knot CI was estimated by assuming the same flux ratio of the knot CI to the knot C at both optical (1.00 x 10^^) 
Hz and X-ray band. 

''The error bars for the PW05 X-ray data were estimated by assuming the same relative precision in both PW05 and M02 data (see text). 




Fig. 1. — Left panel: plot of the SEDs for the knots in the M87 jet from radio through X-ray 
wavelengths. Right panel: plot zooming in on the optical and UV region, where most of 
the data points are and also where most of the curvature is. The data from radio to UV 
arc plotted as filled circles, the X-ray data from PW05 as filled star. The solid lines display 
model fits that include the X-ray data from PW05. The error bars of most measurements 
are too small to be seen here. The best-fit parameters are listed in Table 2. 
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Table 2. Parameters for Model Fits for Radio through X-ray Data 



Parameter 


i/bi(1014Hz) 


z/s2(10i^Hz) 


P 




knot D 


6.21 


97.5 


2.39 


2.34 


knot E 


6.97 




2.40 


3.73 


knot F 


7.30 


10.1 


2.36 


9.33 


knot A 


6.35 


6.94 


2.34 


1.89 


knot B 


4.50 


1.69 


2.35 


4.39 


knot CI 


3.00 


1.34 


2.38 


2.55 



Note. — Col. (1): Knot designation. Col. (2): First 
break frequency in 10^^ Hz. Col. (3): Second break fre- 
quency in 10^^ Hz. Col. (4): Spectral index of electrons. 
Col. (5): Reduced chi square. 



